Multichannel Anomaly of the Resonance Pole Parameters Resolved 
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Inspired by anomalies which the standard scattering matrix pole-extraction procedures have pro- 
duced in a mathematically well defined coupled-channel model, we have developed a new method 
based solely on the assumption of partial-wave analyticity. The new method is simple and appli- 
cable not only to theoretical predictions but to the empirical partial-wave data as well. Since the 
standard pole-extraction procedures turn out to be the lowest-order term of the proposed method 
the anomalies are understood and resolved. 

PACS numbers: 11.55.-m, 11.55.Fv, 14.20.Gk, 25.40.Ny. 



Introduction. The determination of the scattering 
matrix (S-matrix) is considered to be the major objective 
of both, scattering theory and energy-dependent analysis 
of scattering data. The collection of S-matrix poles in the 
"unphysical" Riemann sheet is related to resonance mass 
spectrum 0, Q so obtaining them is the crucial goal of 
any partial- wave analysis. There is, however, a long last- 
ing (and yet unresolved) controversy on the resonances' 
physical properties. It is not clear whether physical mass 
and decay width of a resonance are given by the "conven- 
tional" resonance parameters like Breit-Wigner mass and 
the decay width, or by resonance pole parameters — real 
part and — 2x imaginary part of pole 0,01 . In the case of 
baryon resonances, the compromise is achieved in a way 
that all the conventional, as well as pole parameters are 
collected in the Review of Particle Physics (RPP) 0- 

The motivation for this work came from the fact that 
pole parameters we extracted from our coupled-channel 
partial waves by a standard model-independent pole 
searching method (speed plot 0) were not unique. The 
obtained pole positions varied from one reaction to an- 
other. The time-delay 0, another method of choice for 
resonance pole extraction was employed, but with simi- 
lar outcome. This anomalous behavior challenged com- 
mon sense and the conclusion was drawn that either our 
partial- wave analysis or the applied pole extraction meth- 
ods were incorrect. The extraction methods were care- 
fully examined, and those methods were determined to be 
at fault. This effort resulted in a new model-independent 
extraction method free from this anomaly. In addition, 
we calculated elastic pole residues, and the obtained val- 
ues were in quite good agreement with others published 
in the RPP Since there are no RPP estimates for 
elastic residues yet, this result lends support to provide 
them. 

It is important to draw a distinction between the time- 
delay signal as resonance's "smoking gun" , and utiliza- 
tion of the time-delay (or speed plot) for the extraction 



'Electronic address: lsasa.ceci@irb.hrl 



of resonance parameters. The authors agree that former 
plays a central role as one of resonance existence criteria. 
However, we disagree profoundly with the interpretation 
of the latter as proper resonance parameter extraction 
method. In fact, we show that the speed plot is just the 
lowest-order approximation to a correct pole extraction 
method. 

In this Letter we introduce two methods for obtaining 
the resonance pole parameters from energy-dependent 
partial waves. One is based on analytic continuation of a 
channel propagator. The analytic continuation happens 
to be dependent on the analysis model, since the channel 
propagator is a feature of the Carnegie-Mellon-Berkeley 
(CMB) analysis The other, T-matrix "regulariza- 
tion" (i.e. pole eliminating) method is formulated in the 
model-independent manner, so pole extraction is not re- 
stricted to CMB formalism. Moreover, the only required 
information may be attained from measurable physical 
processes. Therefore we recommend it as The Method 
for obtaining resonance pole parameters. 

The N* pole parameters given in this Letter were 
extracted from partial waves obtained in our current 
partial- wave analysis . 

The Core of CMB Approach. Our current partial- 
wave analysis |^ is based on the CMB approach 
The most prominent property of this approach is ana- 
lyticity of partial waves with respect to Mandelstam s 
variable. In every discussion on partial-wave poles, ana- 
lyticity plays a crucial role since poles are situated in a 
complex plane, away from physical region. Any knowl- 
edge about the nature of partial- wave singularities would 
be impossible to gain if partial waves were not analytic 
functions. The ability to calculate pole positions is not 
just a benefit of the model's analyticity but also a neces- 
sity for the CMB resonance extraction. In this approach, 
the resonance itself is considered to exist if there is an 
associated partial- wave pole in the "unphysical" sheet. 

The central role of the CMB analysis belongs to the 
unitary-normalized partial-wave T-matrix T(z) 0,0,0]. 
It is a matrix in channel indices, and generic com- 
plex variable z in this section denotes Mandelstam s. 
The connection between S- and T-matrix is given by 
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S{z) = I + 2iT{z) where I is the unit matrix. Two 
main ingredients of the model are channel propagator 
$(z) - the diagonal matrix in channel indices which takes 
care of channel related singularities, and bare resonant 
propagator Go{z) - the diagonal matrix in resonant in- 
dices incorporating real first-order poles related to res- 
onances (and background). Background contribution is 
given by two sub-threshold poles (another pole may be 
placed further above considered energy region). Dressed 
channel propagator G{z) is given by the resolvent (Bethc- 
Salpeter) equation G~^{z) = Gq^{z) — S(z), where self- 
energy term S(z) is built from the channel propagator as 
7 • $(z) • 7""". The parameter matrix 7 is a non-square 
matrix obtained from the least-square fit to experimental 
or partial-wave data. In addition to 7 parameter matri- 
ces, the values of the bare propagator real poles arc con- 
currently acquired from the same fit. The partial- wave 
data arc fitted by the unitary-normalized partial-wave 
T-matrix given by relation 



T(z) = v/lm *(z) • 7T . G(z) • 7 • v/lm ^{z). (1) 

Channel propagator matrix $(z) is assembled from 
channel propagator functions (i>(z). The dominant sin- 
gularity in the resonant region, apart from resonances 
themselves, is the physical (channel opening) branching 
point Xo- In the CMB approach, contributions from other 
singularities (left-hand cut, nuclcon pole etc.) are given 
partly by the design of the channel propagator imaginary 
part, while the rest is taken care of by the background. 

Analyticity of channel propagator function 0(z) is en- 
sured by once-subtracted dispersion relation 



Im 0(a;') dx' 
(x' — z){x' — Xo) ' 



(2) 



where P stands for Cauchy principal value. The phys- 
ical (unitarity) branch cut is, thus, chosen to go from 
the branching point Xq to positive infinity. The variable 
x' is used in the integral rather than z' to indicate the 
integration path is on the real axis. 

The form of the channel propagator imaginary part is 
given as 



lm.(j){x) 



[g(x)] 



2L+1 



Qi + \/Ql + [q{^)Y 



2L ■ 



(3) 



where q{x) is the standard two-body center of mass mo- 
mentum for a particular meson-baryon channel, Qi and 
Q2 are model parameters with values equal to the tt me- 
son (or, in our case, the channel meson [6|) mass. L is the 
orbital angular momentum number of the given partial 
wave. 

Extraction Method One: Analytic Continuation. 
From Eq. |j3J) it is evident that 4){z) has a square-root 
type singularity. Instead of calculating the dispersion in- 
tegral jSJ for each point in complex plane, we decided to 



use the expansion (similar to Pietarincn's in Ref. ^3) 

TV 

0i(z) = ^c„(Zi(z))", (4) 

where c„ are coefficients of expansion. The new channel 
dependent variable is given by its principal branch 



Zi{z) 



a - \/xo - z 
a + ^Xo - z ' 



(5) 



with the tuning parameter a. This function is fitted to a 
datasct consisting of imaginary parts of <i){x) from Eq. (O 
and real parts of (\){x) calculated from dispersion relation 
(O, both of them evaluated at real axis (hence x). The 
general idea is that the 0(z) inherits analytic structure 
from Z{z). We obtained parameters a and coefficients 
c„ for each channel, and for all analyzed partial waves. 
The least-square fit is considered to be good if it meets 
following conditions: (i) small number of coefficients c„ 
needed (7 or 8, at most), (ii) the function fitted to the 
part of datasct, when extrapolated outside of the fitted 
region is consistent with the rest of data, and (iii) fitting 
just imaginary part of i>(x) produces real part that is in 
agreement with values obtained from |(2J). 

The channel propagator given by expansion I^J is ob- 
tained quite accurately and works very well in the reso- 
nant region in the vicinity of physical axis. 

Every channel opening is responsible for two Riemann 
sheets: the first (physical) sheet with physical partial 
waves, and the secondary (unphysical) sheet with reso- 
nant poles. To get to the unphysical sheet it is enough 
to use the second branch oi Z{z) 



Zu{z) 



a + ^/Xo ~ z 
a - ^Jxo - z 



(6) 



Finally, it is evident from Eq. that all poles of each 
partial wave must be by construction the same in every 
channel and, in fact, equal to poles of the resolvent G(z). 
Here we use T-matrices obtained in our last partial-wave 
analysis |^ and collect all poles of G(z) obtained by an- 
alytic continuation in Tabled Since this method gives 
poles of partial waves in Mandelstam s variable, compar- 
ison to RPP estimates is made with the square root of 
the Mandelstam pole (selecting branch with positive real 
part) denoted by [i. 

The Anomaly Appears. It is well known that each 
resonance pole on the unphysical sheet is accompanied 
with poles on other Riemann sheets 0, lisl . They are 
attached to the same resonance phenomena, but with 
dissimilar pole positions as well as elastic pole residues 
[Tsl . Pole parameters presented in the RPP Q are taken 
from the unphysical Riemann sheet (the one closest to 
the physical sheet). To be sure that the simple recipe 
given by Eq. © provides us with proper pole parameters 
(i.e. that we are searching for poles on the right sheet), 
we compared results with those obtained by standard 
model-independent pole extraction methods. The two 
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TABLE I: The N* resonance pole parameters obtained by two methods from this Letter along with RPP estimates. The 
N/E is written if a resonance pole position does not have RPP estimate, while the N{1111) stands for resonances unnamed 
in the RPP. For the Af(1535) Sn resonance, the extracted parameters were strongly influenced by the singularity of the rjN 
channel opening so the minimizing process was onerous. 
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Ai'(1535) 


Sii 


1505(10) 


170(80) 


1517 


190 


1522 


146 


19 


-146 


iV(1650) 


Sii 


1660(20) 


160(10) 


1642 


203 


1647 


203 


84 


—58 


Ar(2090) 


Sii 


N/E 


N/E 


1785 


420 










A''(1440) 


Pll 


1365(20) 


210(50) 


1359 


162 


1354 


162 


47 


—95 


Af(1710) 


Pll 


1720(50) 


230(150) 


1728 


138 


1729 


150 


52 


-156 


N{im) 


Pll 


N/E 


N/E 


1708 


174 










Af(2100) 


Pll 


N/E 


N/E 


2113 


345 


2120 


347 


31 


-59 


iV(1720) 


Pl3 


1700(50) 


250(140) 


1686 


235 


1691 


235 


19 


-112 


iV(1520) 


Dl3 


1510(5) 


115(5) 


1505 


123 


1506 


124 


36 


-14 


Af(1700) 


Dl3 


1680(50) 


100(50) 


1805 


130 


1806 


132 


7 


-36 


Ai^(2080) 


Dl3 


N/E 


N/E 


1942 


476 










iV(1675) 


Dl5 


1660(5) 


140(15) 


1657 


134 


1658 


138 


25 


-20 


iV(2200) 


Dl5 


N/E 


N/E 


2133 


437 


2145 


439 


22 


-71 


iV(1680) 


Fl5 


1670(5) 


120(15) 


1664 


134 


1666 


136 


45 


-26 


Af(1990) 


Fl7 


N/E 


N/E 


1990 


303 


2016 


318 


8 


-25 


jV(????) 


Gl7 


N/E 


N/E 


1740 


270 


1749 


280 


6 


-86 


^■(2190) 


Gl7 


2050(100) 


450(100) 


2060 


393 


2068 


389 


34 


-30 



renowned extraction methods are speed-plot 0, [ll [13 
and time-delay 0, . Both methods rely on following 
parameterization of T-matrix elements (Brcit-Wigner) 

nz)^ ^ + (t{z)-^\, (7) 
fi — z \ fx — z J 

resonant part smooth background 

where ji and r are pole position and pole residue, respec- 
tively. Here z stands for center-of-mass energy (v^). If 
one plots modulus of the "speed" of T (i.e. \dT [z) / dz\) , 
the resonance produces a peak in this speed plot. There 
are, however, known exceptions like A^(1535) which is 
hidden "under the cloak" of the rjN channel opening Q . 
Resonance poles are extracted from Eq. Q under the 
erroneous assumption that "speed" of the background 
can be completely neglected when compared to "speed" 
of the resonant part. We compared the speed-plot pole 
positions of TriV elastic T-matrix element to positions ob- 
tained by using other T-matrix elements: from other 
quasi-elastic (like i-jN rjN) as well as inelastic pro- 
cesses (e.g. -kN ilN). Obtained pole positions were 
shifted by a few tens of MeVs which indicated that there 
was something wrong. 

The time-delay is given as the imaginary part 

of the product of S-matrix inverse and S-matrix "speed" 
(Im [5'l^(z) • dS{z)/dz'\). The resonance reveals itself as 
peak in a plot of time-delay versus z. The parameters are 
again obtained by neglecting background contributions 
thus discrepancies that occurred were similar (though 
less) to those of speed plot. 

The background contribution can be disregarded near 



a pole. However, the real axis may be too far from 
the pole, so observed discrepancies can be explained by 
background contribution. Instead of carefully designing 
this background, we developed a following method to ap- 
proach the complex pole (where background is negligible) 
without leaving the real axis. 

Extraction Method Two: Regularization. Let there be 
an analytic function T{z) of complex variable z that has 
a first-order pole at the some complex point /i. This 
function T{z) can be any T-matrix element, and variable 
z can be either Mandelstam s or center-of-mass energy 
^/s. The latter is used in this method in order to achieve 
full correspondence with the speed plot. Since all physi- 
cal processes occur at real energy values, we are allowed 
to determine directly only values of T{x), where x is a 
real number. To be able to successfully continue our T{x) 
into the complex plane to search for its pole(s), we should 
"regularize" this function (i.e. somehow remove the sin- 
gularity). Then, any simple expansion would converge in 
the proximity of the removed pole. The appropriate way 
to "regularize" the function with a simple pole at fj,, is to 
multiply it by a form that has a simple zero at the same 
point 

/(z) = (M-2)T(z). (8) 

From this definition and Eq. {Tj), it is evident that the 
value of f{fj,) gives the residue r of T{z) at point fi. As 
we have access to the function values on real axis only, 
the Taylor expansion of / is done over some real x to 
give the value (residue) in the pole fi (where background 
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is highly suppressed) 



N 



x)" + Rn{x, m)- 



(9) 



n=0 



The expansion is exphcitly written to the order N, and 
the rest is designated by i?Af (a;, /^). The nth derivative of 
f{x) in the form of T(x) is given as 



/(")(a;) = (/i-a;)r(")(a;)-nr("-i)(x). 



(10) 



Insertion of this derivative into Taylor expansion conve- 
niently cancels all consecutive terms in the sum, except 
the last one 

f{^^) = ^^^^{^^-xy''+''>+R^{x,^,), (ii) 

where T^^^{x) is the A^th energy derivative of T-matrix 
element. To simplify the notation, the pole can be writ- 
ten as some general complex number jj, = a + ih. Once 
Taylor series converges (i.e. the rest Rn{x,h) is disre- 
garded) the absolute value of both sides of Eq. Hll|l yields 
to 



I/(m)I 



|rW(a;)| 

m 



\a + ih — x\ 



(TV+l) 



(12) 



To keep the form as simple as possible Eq. (|12|l is raised 
to the power of 2/{N +1). The elemental second-order 
polynomial emerges after simple rearrangement on one 
side of the equation 



{a - xf + &2 



i/(/^)r 



(TV! 



(13) 



where the part that is attainable from energy-dependent 
analysis is put on the right-hand side. The fitting func- 
tion with just three fit parameters a, h and |/(Ai)| is writ- 
ten on the left-hand side. The speed plot method turns 
out to be identical to the first order of this relation (with 
N = 1). 

The dataset was produced by the right-hand side of 
Eq. (HSJ). rW(a;) was given by numerical derivation of 



energy-dependent partial waves obtained in our analysis 
@. A step of 2 MeV provided a stable procedure. The 
Taylor expansion is considered to converge when the ex- 
tracted parameters settle down. The higher orders were 
used to obtain more accurate values of pole parameters. 
To acquire reliable fit results, we considered data grouped 
in a parabolic shape (in accordance with the second-order 
polynomial) . 

The elastic-pole residue is commonly given by its 
absolute value \r\ and its phase 9, namely 

|r| = |/(M)|, tan0 = Im/(Ai)/Re/(/i), (14) 

where this particular (elastic) /(/i) is given by the first 
term in Eq. Hll|l with T{z) being nN elastic T-matrix 
element. 

Pole parameters attained in this way from -kN elastic 
process are given in Table ^ In order to verify the proce- 
dure, we applied it to other channel processes. Contrary 
to anomalous results obtained when using standard pro- 
cedures, inelastic poles varied by only a few MeV from 
the elastic ones. 

Conclusions. Our new analytic continuation method 
(Pietarincn expansion of the channel propagator) pro- 
vides pole positions quickly and precisely while avoiding 
problems with numerical principal value integration and 
interpolation. The obtained pole positions arc in accor- 
dance with RPP values. 

The detected anomaly of standard speed plot and time 
delay methods for pole extraction is fundamental because 
the unknown contribution (energy derivative of back- 
ground on the real axis) is considered to be insignificant. 
The new "regularization" method successfully finds res- 
onance pole parameters from a T-matrix in a model in- 
dependent way. The advantage of the given method is 
that it can generally be applied to most analytic func- 
tions that have a simple pole and values known on any 
line segment reasonably close to the pole. 

The elastic pole residues arc in accordance with those 
given in RPP. Since there are still no residue estimates in 
RPP, we strongly advocate making them for the future 
editions. 
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